Dissipation in ferrofluids: Mesoscopic versus hydrodynamic theory 
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Part of the field dependent dissipation in ferrofluids occurs due to the rotational motion of 
the ferromagnetic grains relative to the viscous flow of the carrier fluid. The classical theoretical 
description due to Shliomis uses a mesoscopic treatment of the particle motion to derive a relaxation 
equation for the non-equilibrium part of the magnetization. Complementary, the hydrodynamic 
approach of Liu involves only macroscopic quantities and results in dissipative Maxwell equations 
for the magnetic fields in the ferrofluid. Different stress tensors and constitutive equations lead to 
deviating theoretical predictions in those situations, where the magnetic relaxation processes cannot 
be considered instantaneous on the hydrodynamic time scale. We quantify these differences for two 
. . - situations of experimental relevance namely a resting fluid in an oscillating oblique field and the 

damping of parametrically excited surface waves. The possibilities of an experimental differentiation 
between the two theoretical approaches is discussed. 

PACS: 75.50.Mm, 47.32.-y, 47.35.+i 

D ' 47.35.+i Hydrodynamic waves 
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The interplay between hydrodynamic behavior and magnetic field sensitivity gives rise to a variety of new physical 
effects and makes ferrofluids a fascinating subject of research with many interesting technical applications [Q-^. 
_ Practically all hydrodynamic properties can be modified or modulated by an external magnetic field. Although very 
0^ ■ similar effects are possible with electrically polarized fluids, the coupling of the fields to the hydrodynamics is more 
[ pronounced in the magnetic case. Since the experimental handling of strong magnetic fields is simpler than that of 
electric ones, ferrofiuids are the experimenters choice to analyze the field affected hydrodynamic motion. 

On the other hand a thorough theoretical understanding of the many aspects of ferrofluid motion is a rather 
Q . ambitious task |^,^. This is in particular true if the motion is such that dissipation effects are important. The 
c/3 ' theoretical analysis must then account for the deviations from and the relaxation towards equilibrium in a consistent 
fjj . and quantitatively accurate way. 

One of the interesting properties of ferroffuids is their tunable viscosity which can be controlled by an external 
O |, magnetic field. This so called magneto-viscous effect has first been observed experimentally by Mc Tague Quali- 
tatively it can be accounted for 0| by magnetic torques acting upon the suspended ferromagnetic particles: depending 
on the relative orientation between magnetic field and local vorticity the particle rotation is hindered. The friction 
at the coated particle surfaces generates an extra dissipation and thus leads to an enhanced effective viscosity. This 
^ . additional rotational friction is described by a "rotational viscosity" rjji. Based on the theory of rotational Brownian 
■ " ' motion of non- interacting rigid dipoles, Shliomis developed a theoretical description of this effect. The basic ingre- 
dient is a relaxation equation for the non-equilibrium part of the magnetization resulting from a stochastic description 
of an ensemble of ferromagnetic grains. This treatment allows to derive an analytic expression for the magnetic field 
dependence of rjji 0|. We will refer to this approach as the mesoscopic theory since it starts with a rather detailed 
characterization of the processes occurring on a 10 nm scale in order to determine macroscopic properties like the 
viscosity. 

Quite complementary the same phenomena can be described within the framework of non-equilibrium thermo- 
dynamics. Appropriate magnetic field variables must then be introduced into the set of relevant thermodynamic 
quantities. New dissipative currents show up when the system is driven out of equilibrium. They are coupled to 
the thermodynamic forces by additional Onsager coefficients. This program was carried out by Liu and coworkers 
in a series of papers 1^,^,^ and gave rise to what is called the hydrodynamic Maxwell theory. Making no reference 
at all to the microscopic mechanism of dissipation this approach has the appealing feature of being very general and 
applicable to electrically or magnetically polarizable continuous media. As in every hydrodynamic theory the Onsager 
coefficients are free parameters which have to be determined experimentally or from a microscopic theory. 

The detailed quantitative relation between the two theoretical approaches is presently not completely clear and 
there has been some controversy in the literature over the recent years. In order to contribute to a clarification of 
this issue we investigate in the present paper the consequences of both the mesoscopic and the hydrodynamic theory 
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for two simple experimental setups involving ferrofluids. On the one hand we analyze the influence of an oscillating 
oblique magnetic field on a ferrofluid at rest, on the other hand we investigate the damping of parametrically driven 
surface waves on the ferrofluid in a constant external magnetic field perpendicular to the undisturbed surface. Both 
situation are easily accessible to experiments. Interestingly the two theoretical approaches give rise to different results 
for some of the relevant quantities. We finally discuss whether these differences are pronounced enough to allow an 
experimental decision between the theories. 



II. EQUATIONS OF MOTION 

The fluid motion is governed by the hydrodynamic balance equations for mass and linear momentum. For an 
incompressible ferrofluid (density p) with a velocity field v(r,t) these equations read as 

V-v =0 (1) 
p{dt^r + V V)v = V ■ n + pg, (2) 

where 11^ is the stress tensor containing reactive and dissipative contributions (see later) and g = (0, 0, —g) is 
the acceleration due to gravity. The two theoretical approaches under consideration differ in their expressions for 
the stress tensor and the field equations. It will turn out that they produce different results when the magnetic 
relaxation time r (usually between 10~^ and lO^^s) cannot be neglected with respect to the time scale, l/w say, 
of the hydrodynamic motion. In those cases it is argued that the local magnetization M(r,t) deviates from its 
local equilibrium value Meg(r,t) giving rise to an additional dissipation, which we shall denote here as "magneto- 
dissipation" . This phenomenon is responsible for the field dependent viscosity , the magneto- vortical resonance 
1^,^, or the spectacular negative viscosity effect [^-^ detected recently in ferrofluids. 

In order to simplify the mathematics we make throughout this paper the following approximations which do not 
remove the basic differences between the two approaches: (i) The magnetic field H within the ferrofluid is weak 
enough to warrant a linear relationship Mg^ — xH, where the susceptibility is supposed to be proportional to the 
density p. (ii) The local non-equilibrium magnetization M deviates only slightly from Meg. (iii) Magneto-dissipation 
is treated in the so called low-frequency limit characterized by tlo <^ 1. 



A. Field equations in the mesoscopic theory 

Following Shliomis 0| (see also Ref. the stress tensor for ferrofluids appears in the form 



Trmes 
ij 



H^)S,j+H,Bj+r]{y,Vj 



Mo 



e,,fe(MxH)A 



(3) 



In Eq. (||) occurs the pressure field p{r,t), the Maxwell stress tensor with the induction B = /io(H + M) and the 
usual dissipative contribution proportional to the shear viscosity rj. Moreover, the last term on the right hand side 
of Eq. (^ arises if the magnetization exhibits an off-equilibrium component perpendicular to H. It is related to 
the magneto-dissipation and guarantees the symmetry of the stress tensor if the local directions of H and B do not 
coincide. The appearance of M in the stress tensor (|3|) requires an extra constitutive equation for the off-equilibrium 
component (5M = M — xH of the magnetization. Under the above simplifications (i-iii) the following relation can be 
established [01 



SM ^ -rx 



dtU + (v • V)H - ZilZ X H 



where r denotes a relaxation time. 

The magnetic part of the problem is treated in the framework of the magnetostatic approximation 

V X H = 
V • B = 0. 



(4) 



(5) 
(6) 



With this simplification and the above expression for the stress tensor (||) the hydrodynamics obeys the Navier-Stokes 
equation, which assumes the form 



Mo, 



pdtv + p(v • V)v = -Vp + pg + T]V V + po{M. ■ V)H + ^^V x (M x H) 



(7) 
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B. The field equations of the hydrodynamic theory 

The complete set of hydrodynamic Maxwell equations is given by 

9tD = VxH~je/ a^B^-VxE (8) 
V ■ D = V • B = 0, 

where E and D denote the electric and displacement fields. For non-conducting ferrofluids there are no electric charges 
or currents, thus je/ = and pei =0. The fields H and E split into reactive and dissipative parts 

H = + and E = E^ + E^ 

where = de/ffR and E^ = de/dD are derivatives of the thermodynamic energy density e and contain only 
equilibrium information. The dissipative fields and E^ represent the off-equilibrium currents and are functions 
of the thermodynamic forces. For isotropic isothermal situations and neglecting off-diagonal couplings they depend 
only on the electromagnetic fields in the local rest frame and are given by |^ 

H^ = _iiVx (E^ + vxB) and E^ - V x (H^ - v x D), (9) 
Mo eo 

where the parameters a and (3 are Onsager coefficients. In the present paper we restrict ourselves to ferrofluids 
with a negligible dielectric susceptibility Xei — (e.g. hydrocarbon based ferrofluids) giving rise to /3 ~ (see also 
[^) implying E^ — 0. Moreover, the ratio between time and space derivatives of fields respectively is uj/(ck) and 
for hydrodynamic frequencies uu ^ lO^s^^ and wavelengths 1/fc ~ lQ~^m~^ we may neglect the time derivative of 
D in the Maxwell equations (||). This means that the magnetic field may be assumed to follow the fluid motion 
instantaneously. In this way we recover the magnetostatic field equations 

V-B =0 (10) 
V x H = . 

At a liquid-air interface these fields are subjected to the usual continuity conditions for the normal component of B 
and the tangential component of H. 

Using incompressibility we get from (^) and (||) the following constitutive relation for the dissipative field 

= — [^tB + (v • V)B - (B • V)v] , (11) 
Mo 

which is the analogue to (5M in the approach of Shliomis. For later use we introduce the reactive contribution to the 
magnetization 

= — B - . (12) 

Mo 

This definition differs from the "real" magnetization M = B//io — H in using only the reactive part H-'^ of the 
magnetic field. The advantage of this formulation is that for the present case of constant and scalar susceptibility the 
fields B,H^ and = are all parallel to each other |jl^. The fact that H and hence M are not parallel to B 

is the very reason for the magneto-dissipative effects [^HJ^ investigated here. 

With the assumption of the susceptibility being proportional to the density p, the stress tensor of the fluid simplifies 
to |,0 



Y^yd _ 



p+>^{H'^f\8,j+H^B,+i^{V,v, + V,Vj) . (13) 



The fluid dynamics is therefore governed by the Navier-Stokes equation, which using ( p^ acquires the form 

pdtv + p(vV)v = -Vp + pg + riV^v + noM^VH^ + B x (V x H^). (14) 

The fourth term on the r.h.s. of the Navier-Stokes equation was simplified to involve only the magnitudes of the 
vectors and by exploiting the fact that they are parallel. 
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C. Comparison between the two approaches 

Let us first emphasize tliat tlie two descriptions introduced above are fully equivalent if the magnetic relaxation 
processes can be considered to be instantaneous (t = a = 0) on the time scale of the hydrodynamic motion. In this 
case of vanishing magneto-dissipation we have (5M = = 0, i.e. M-'^ = M || H = H^, and thus the stress tensors 
and (|l^) or respectively the Navier-Stokes equations and ( p^ coincide. 

Differences between the two approaches arise due to the treatment of the magneto-dissipation occurring at finite 
(5M and H^, respectively. By comparing Eqs. (|^) and (^3|) and expressing the antisymmetric part of the tensor 
element HfBj in @ by (l/2)eyfe(H^ x B)^ = -(^o/2)eijfe(M x U)k one arrives at 

^hyd _ jjmes ^ ^^jj , _ \{H^ B, + HfB,), (15) 

where terms of order [H^Y have been dropped. The stress tensors in Eq. ([l^) differ by a diagonal element and a 
symmetric non-diagonal contribution. The isotropic term just leads to a re-normalization of the pressure, which we 
won't pursue here. The non-diagonal tensor, however, gives rise to a dissipative force which might be detectable in 
an experiment (see below). 

It is also instructive to compare the constitutive equations, which determine the dissipative fields: From the identity 

/^o(H + M) = B = /io(H^ + M^) = Ato(l + x)H^ (16) 
one finds that the dissipative field contributions (5M and coincide up to a pre-factor: 

H^ = 5— (M-xH) = ^—Sm. (17) 



The associated constitutive relations 

5M. = -Tx 



9tH + (v ■ V)H - i(V X v) X H 



(18) 



= ii[atB + (v V)B-(B- V)v], (19) 
Mo 

however, are not fully compatible: Omitting contributions of 0(q;cij)^ Eq. (19) can be re-casted as 

= a(l + x)[dM + (v • V)H - (H ■ V)v]. (20) 

By comparing the two first terms on the right hand side of Eqs. (18, |20| ) and using ( |l7| ) one can relate the characteristic 
relaxation times by ||l^ 

- - ^ (21) 



(i + x)2- 

The third terms in Eqs. (18j20|) describe different physics. This is easy to see by considering a homogeneous stationary 
magnetic field: According to the mesoscopic theory, Eq. (18), only a rotational flow is able to drive a magneto- 
dissipative field. Contrary in the hydrodynamic approach, Eq. (|2^) states that a flow gradient parallel to the magnetic 
field suffices to excite a finite . It is interesting to point out that the deviation between the constitutive equations 
Eqs. (18, ^0|) drops out for a solid body rotation flow field v = u; x r, which is often investigated in the ferrofluid 
literature. 

To summarize this section: We observe two differences between the classical description of Shliomis and the hy- 
drodynamic Maxwell theory of Liu: One in the stress tensors Ily and a second in the constitutive equations for 
the dissipative fields (5M and H'^ respectively. In the following we investigate the possibilities of an experimental 
differentiation between the two theoretical approaches. 



III. FERROFLUID AT REST IN AN OBLIQUE OSCILLATING MAGNETIC FIELD 

We consider a resting ferrofluid layer exposed to a linearly polarized magnetic field H^^' — {H^\q^ h'^'^) con ojt 
directed oblique to the surface. Here and in the following flelds in the air carry the superscript " (2)" , while fields 
within the ferrofluid are denoted without superscript. 
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The stress tensor in the air above the ferrofluid is taken as 



aim _|_ A^O ^jj(2)-^2 



(22) 



where we assume a constant atmospheric pressure p"*"*. 

For a fluid at rest the constitutive equations (18,19) coincide and thus eventual discrepancies between the two 
theories wiU result from the deviating stress tensors. To test possible differences we evaluate the tangential force 



Ilxz acting upon a free surface element. By using the magnetic interface conditions for H and B and 



owing to the fact that isotropic tensor elements do not contribute to tangential forces we obtain 



and 



Mo 



(24) 



(25) 



The computation of ft requires knowledge of H and B within the ferrofluid. From the respective constitutive equations 
and the interface conditions we flnd 







H 



(2) 



1 + X - ii^TX 



e'"* + c.c. 



■Qhyd ^ Mo 

2 



H 



(2) 



1 

L i+x 



iLua 







Evaluating the tangential stresses ( p^ , p5| ) up to leading order in tuj or respectively auj we get 

= O[{TL0f] 



and 



rhyd Mo 
ft = Y (1 



x)hPhP sm2ujt + 0[{aujf 



(27) 
(28) 

(30) 
(31) 



At the considered accuracy level, the mesoscopic approach states that the resting ferrofluid is in equilibrium, while 
the hydrodynamic Maxwell theory predicts a residual tangential surface force oscillating with twice the excitation 
frequency. This force drops out in the static limit ^ or when the applied field is directed normally or tangentially 
to the surface. The maximum effect is achieved for an inflection angle of 45°, where the product of the incident held 

(2) (2) 

components Hx Hz is maximum. We conclude that the resting fluid cannot be in equilibrium, i.e. the force will 
drive a convective motion. With the aid of appropriate tracer particles exposed to the surface of the fluid we expect 
the acceleration to be easily detectable in an experiment. 



IV. SURFACE WAVES WITH MAGNETO-DISSIPATION 

In this section we investigate a setup, which aims at probing the different constitutive relations (18,19). As outlined 
in Sec. II C a non-rotational flow profile is most appropriate for this purpose. Surface waves on a deep fluid are a 
classical example: The associated flow field is purely potential, except in a thin boundary layer along the system 
boundaries. By increasing the system dimensions this parasitic influence can be limited to the viscous skin layer 
beneath the free surface. Our aim is to derive the complex dispersion relation for low frequency small amplitude 
surface waves with magneto-dissipation taken into account. 



A. Surface waves without magneto-dissipation 



In this subsection we review the present knowledge on small amplitude surface waves on a ferrofluid in a static 
magnetic field Hext normal to the surface. With z — Cix,t) oc e^i'^^-'^*-) we denote the surface deflection in vertical 
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z-direction depending on the horizontal coordinate x and on time t. The real part of the complex valued function 
uj{k) = uj'{k) + iuj" [k) reflects the wave dispersion, while the imaginary part accounts for the decay rate. 

For a non- viscous [rj = 0) ferrofluid Lij"{k) vanishes identically and the inviscid wave dispersion (indicated by the 
subscript " 0" ) is given by (l9[0] 

a;2(fc, H,,,) ^gk- ^ . Hi, k' + lk\ (32) 

P (1 + X)(2 + X) P 

Here 7 denotes the coefficient of surface tension. The magnetic field leads to a negative contribution proportional to 
k^. If Hext is sufficiently strong luq becomes negative, indicating the onset of the Rosensweig instability. 

The effect of viscosity on the dispersion has been investigated theoretically in [^,0 under the assumption of an 
infinitely fast magnetic relaxation (vanishing magneto-dissipative effect). It is found that the expression for u}{k) can 
no longer be given explicitly, rather the following implicit relation applies, 

V{k,uj) =uj^ -Lo^{k)+X{k,uj) = (33) 

where the viscous contribution reads 

X(fc,w) = 4iwiyfc2 + 4(iyfc2)2{i _ 1}. (34) 

k 



Here v — rj/ p is the kinematic viscosity and q = ^Jk"^ — iuj /v. Note that the viscous contribution does not depend 
on the magnetic field and therefore coincides with the expression for non- magnetic fluids p3]. This is because the 
magnetic contributions to the stress tensor are purely conservative at vanishing or H . In the mathematical 
derivation of (33 Sj) this is reflected by a decoupling of the magnetostatic from the hydrodynamic problem. Later we 



will see, that this simplification no longer applies if magneto-dissipative corrections are taken into account. 



B. Formulation of the problem with magneto-dissipation included 

The investigation requires solving the equations of motion together with appropriate boundary conditions at the free 
liquid air interface. We consider a horizontally unbounded ferrofluid layer of inflnite depth (from z = to z = —00) 
in contact with air at its free surface. The layer is exposed to an external stationary magnetic field "H-ext — (0, 0, Hext) 
directed perpendicularly to the undisturbed surface [see Fig. 1]. 

The balances of the normal and tangential stresses at the liquid air interface are expressed by the conditions 

n,(Hij - Iifj')nj = -7(Vjnj)n, (35) 
t,(H,, -ng))n, =0 (36) 

where t(r, t) and n(r, t) denote local unit vectors tangential and normal to the surface respectively. The stress tensor 
in the air above the ferrofluid is given by Eq.(|2^). The discontinuity of the normal stress condition arises by virtue of 
the flnite surface tension 7. The kinematic surface condition requires 

dtC, + udxC = w, (37) 

where v = {u, 0, w). The second lateral space direction y need not to be considered. 
As usual the boundary conditions for the magnetic field at the free surface read as 



(B-B(2)).n = 0, 
(H - H(2)) . t = 0. 



(38) 
(39) 
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C. The motionless state and linearized equations for the perturbations 

As reference state (index 0) we consider the motionless (v — 0) fluid with a flat (C = 0) surface. The associated 
magnetic quantities are 

Bo =B^') = (0,0,Aioife.t) (40) 

Mo = (0,0, -^i7e.t) 
l + X 

l + X 

H^'^ = (0,0,i7e.t) 

= (5M = 0. 

^From the momentum balance with the appropriate boundary condition we get for the pressure 

2 



Po{x,z} = -pgz- — 



X 



l + X 



HLf (41) 



Next we study the time evolution of small perturbations by linearizing the field equations and their boundary condi- 
tions around the basic state solution. To this end we introduce small deviations denoted by v = (u,0, w), 5p^ ^, b, 
h^'^\ h, h'^^\ m. Linearizing the Navier-Stokes equations (0) and ([l^ ) we respectively find for the mesoscopic and the 
hydrodynamic approach 

pdtv = -V6p + riAv + fioMod^h + y V x [Mq x h + m x Hq] (42) 

and 

pdtv = -\/Sp + 77AV + MoMoV/i, + /ioi?erf(VHf - 9,H^). (43) 

Here we have exploited the fact that Hq and Mo are both parallel to the z-dircction. The linear magnetostatic field 
equations remain unchanged 

=V-b = V-b(2) (44) 

= V X h = V X h(2). (45) 

The boundary conditions ( ^ , ^9|) are linearized by using the lowest order expressions 

dC dC 
n=(-^,0,l) and t = (l,0,^) (46) 

to get 

6.^5(2), (47) 

h,. - = Mo|^ (48) 
for the magnetic interface conditions. The kinematic surface condition simplifies to 

dtC = w, (49) 

while we get 

= -pgC + Sp + fioMom^ ~ 27jd^w + -fd^C (50) 
for the normal stress balance at the surface. Introducing the dimensionless magneto-dissipative perturbation parameter 

— -l^ext — ~~r^ ; ^^ext- (51) 
77 ?7(l + x) 
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the condition for the tangential stress reads in the mesoscopic approach 

^ {dxw - dzu) = (a^u + dxw) (52) 

whereas for the hydrodynamic Maxwell theory we simply find 

BoH^ =v{dzu + d^w). (53) 

To elucidate the time evolution of the perturbed state we assume for all lateral field dependencies a plane wave 
behavior oc exp[i{kx — cot)], e.g., v = v(z) exp[i{kx — cut)] or ^ = (^exp[i{kx — i^t)], etc. The associated z-dependence 
(if appropriate) is indicated by a bar. We now first determine the magnetic fields from (^,^) in terms of the surface 
deflection C, and the flow field v using the boundary conditions (^ , ^ ) . Then we solve for the velocity field and impose 
the hydrodynamic boundary conditions (^) and (|53|), respectively. The condition for non-trivial solutions of (|50| ) 
finally yields the desired complex dispersion relation a;(fc). The explicit calculations are sketched in the appendix. 

D. Results 

The expressions we find for the implicit dispersion relation T>(k,uj) are of the form 

= 2?'"^"(t^, k)=u? - uolik) + Uujvk^ + 4(zyfc2)2{l - 1} - ti(y}e- f^^{l - 1} (54) 

k 2 + X k 

= V'^y'i.., k) {u? + 2zc.^fc^) ^~^^[;"f ^ - ^l{k) + 2(2 + «)^2^3^1±^ (55) 



lUJ lUJ 



with 



f = fc2 - , , , (56) 

We analyze these results in the following for a ferrofiuid with a magnetic susceptibility and viscosity appropriate for 
surface wave experiments (e.g. APG J12 j23j, Ferrofiuidics) . In a moderate external field up to Hg^t = 15k A/m, 
which is below the Rosensweig threshold, the constant susceptibility approximation holds with an error better than 
4%. At this field strength the perturbation parameter given by Eq. (|l|) is smaller than 0.06. For a quantitative 
investigation of the magneto-dissipative effect it is therefore sufficient to expand I?(/c, w) up to first order in k giving 

V{k,u;) =uj^ -ujl+X{k,u;) + KXM{k,uj,Hext) + 0{K^) (58) 

with the magneto-dissipative contributions 

Xir(fc,^,i^e.) = (^fcY{^-^{f-l}} (59) 



in the mesoscopic theory and 



XJ {k,uj,Hext) = -JTY ^ ' k 



X ■ ^k 



2i- 



(60) 



in the hydrodynamic approach. We observe a different scaling behavior with the shear visosity v. The magneto- 
dissipative correction according to Eq. ( |5S| ) scales with i/"^/^, while it is only proportional to v in Eq. (^0|). This 
deviating proportionality traces back to the constitutive equations. According to Eq. (18) magneto-dissipation is 
associated with rotational flow and thus confined to the thin boundary layer beneath the surface, where damping is 
proportional to v'^/'^ [Q. This is in contrast to Eq. (19), where the dissipative field is finite over the whole fluid 
layer. In particular, dissipation within the convective bulk, which is known [p5|| to scale with provides the leading 
contribution to Eq. (|60|). 
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Due to the smallness of k the magneto-dissipative contribution KXM{k,uj) modifies the ordinary viscous shear 
damping X{k,uj) only sHghtly. Moreover, when studying the field dependence of the damping rate to" {k, Hext) one 
has to account for appreciable wavenumber shifts resulting from the non-dissipative magnetic contribution 



H'Ltk'. (61) 



Mo X tt2 ,2 



p (l + x)(2 + x) 

in 'D(k,Lo). Since X(k,Lo) in turn is largely fc-dependent, the reactive term ( |6l|) strongly feeds back into the effective 
field dependence of u;"{Hext, k) and therefore masks the tiny magneto- viscous contributiony Xm- To circumvent this 
difficulty one has to adapt the wavenumber k such that WQ(fc, Hext) as given by Eq.(^) is held constant when varying 
Hext- This can be accomplished by simultaneously tuning the wave excitation frequency to' during the Hextscan. 
Fig. 2a depicts the field dependence of uj" evaluated according to this protocol. As expected, the damping rate 
increases with the field strength. The predictions of the two approaches deviate by less than 4% (Fig. 2a, b). With 
regard to the approximations made and the uncertainty of the material specifications this tiny difference is presumably 
too small to be resolved in a surface wave experiment. For completeness we show in Fig. 2c, d how the excitation 
frequency uj' must be adapted to guarantee a constant loq during the variation of the field intensity. 



V. SUMMARY AND DISCUSSION 



The aim of the present paper was to compare two different theoretical approaches to the problem of magnetic field 
dependent dissipation in ferroffuids. On the one hand we considered the theory given by Shliomis resting on a meso- 
scopic treatment of the rotation of the magnetic particles. On the other hand we applied the hydrodynamic Maxwell 
theory advocated by Liu which treats field dissipation within the standard framework of irreversible thermodynamics. 
For experimental arrangements in which the magnetic relaxation processes can be considered infinitely fast the two 
descriptions are equivalent. This is, e.g., the case if the ferrofiuid is composed of ferromagnetically soft particles, 
where the associated magnetic relaxation is usually very rapid as it is dominated by the fast Neel process. 

On the other hand, for ferrofiuids composed of materials with a high specific magnetic anisotropy (e.g. elementary 
cobalt or cobalt-ferrites) the magnetic relaxation process is largely dominated by the Brownian rotation of the particles 
in the carrier matrix. In high viscosity carrier fluids the associated Brownian relaxation time can be large and reach 
the time scale of the hydrodynamic motion. This is the situation in which magneto-dissipation must be accounted 
for. We have shown that this effect is treated differently by the approaches of Shliomis and Liu due to the use of 
different stress tensors as well as distinct constitutive equations for the dissipative fields. 

In a ferrofluid layer at rest the constitutive equations coincide. If the surface is exposed to an oblique linearly 
polarized magnetic field, the differing stress tensors lead to different statements about whether the resting ferrofluid 
is in equilibrium or not: Whereas the mesoscopic approach gives no deviation from equilibrium the hydrodynamic 
theory predicts an onset of convection. We expect these predictions to be easy to verify or falsify in an experiment. 

On the other hand, an experimental check of the validity of the constitutive equations is probably more subtle. The 
magneto-dissipative contribution is usually very small as compared with the ordinary shear dissipation. Clearly, it 
becomes most pronounced if the ferrofluid is in solid body rotation, where shear dissipation drops out. Unfortunately, 
this special flow geometry - when introduced into the constitutive equations of Shliomis and Liu - leads to the same 
dissipative flelds. We therefore looked for a non-rotational flow geometry, which is easy to realize and which fulfills 
the requirements of small amplitude low frequency motion. As a possible candidate we propose the flow field of free 
surface waves. On the basis of the two theoretical approaches at hand we worked out the associated complex dispersion 
relations. In the final expressions for the magneto-dissipative contribution we observe a different scaling behavior with 
the shear viscosity v, which could be traced back to the deviation between the constitutive equations. On the basis 
of the material specifications for a real ferroffuid and the experimental parameters for a low frequency surface wave 
experiment we computed the damping rate and evaluated the difference between the predictions. Unfortunately, the 
resulting deviation is less than 4% and therefore smaller than the expected experimental error. 
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VI. APPENDIX 

In this appendix we sketch the derivation of the dispersion relations for surface waves within the two theoretical 
approaches. 



A. The mesoscopic approach 



Due to ( p^ we introduce the scalar magnetic potentials by 

H ^ -Vi/', H(2) 

which obey the following relations 



(2) 



XT 



2(1 + X) 



HqV^w 



vV^'^ = 0. 



(62) 



(63) 



Eq.(|63|) has been derived by using (^ and the constitutive equation (18). Invoking once more the low frequency 
approximation, the latter occurs in the form 



m — xh + (5M = X 
The potentials and are solved by 



h+ -(V X v) X Ho 



(2) 



2(1 + X) 



(2) —kz i{kx—ujt) 



The boundary conditions ( ^ , ^8|) yield the as yet unknown pre-factors 



X 



2 + X 
1 + X 



-Afn 



Ha 



■X 



Mo 



2 + X 

——Hovk (1 - -) 

2 + X K 



vk^{l — —) — iuj 



c, 



(64) 

(65) 
(66) 

(67) 
(68) 



where q ^ -^/fc^ — ioj/[i/(^l + k/4)]. ^From Eq. (65) and the expressions for the coefficients (67,68) it can be seen 
that the solution of the magnetic field equations couple to the velocity field w. This coupling is removed in the 
absence of the magneto-dissipation (r = 0) and the result of the standard inviscid calculation [Q is recovered. These 
considerations complete the solution of the magnetic field problem in terms of the hydrodynamic variables. 

We now turn to the solution of the hydrodynamic problem. ^From the linearized equations 43) we eliminate 
the gradient terms by operating twice with curl on the Navier-Stokes-equation. Projecting on the z-axis and using 
again Eq. (64) we find 



pdt - (1 + ^)?7V2j V^w ^ 0. (69) 
Introducing again the plane wave ansatz for the lateral dependence of w, the vertical dependence is determined by 

{dl-e){d^,-f)w^O (70) 

with 



The proper solution of (|7^) remaining bounded for z ^ — oo is 

w{z) = Ae''' + Be^' 



(71) 



(72) 



with the so far undetermined coefficients A and B. It is always understood that the real part of q is positive. Using 
the balance equation for normal stresses (^0|) and the kinematic surface condition (^9|) allows for to express A and B 
in terms of resulting finally in Eq. (54). 
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B. The hydrodynamic approach 



Again (^5|) allows to introduce scalar magnetic potentials 



H(2) = -V^/>(2) , (74) 



The linearization of Eq. (O) is in the present case 



= — 5tB - aiJext^zV . (75) 
Mo 



^From this equation and (44) we find that due to the incompressibility of the flow the magnetic potentials must be 
harmonic. Therefore 

V> = V; e''^ e'(*="=-"*) (76) 

^{2) ^^i2)^-kz^^{ka:-ut)^ (77) 



Using the boundary conditions (47 42) we find 

_ jHrHextdzw{0) - (1 - iiirauj)C,Mo 

tp = — : 78 

\ + — IfirCtU! 

-(2) _ jfJ,rHextdzW{0) + ^rCMp 

V -I . ■ ) \''^) 

1 + /i,- — tfirauj 

where ^^^(O) denotes the derivative of w{z) at z = and — ^ + X- Again it can be seen that the solution of 
the magnetic field equations couples to the velocity field w. In the absence of the magneto-dissipation {a = 0) this 
coupling is removed and the result of the standard inviscid calculation ||l[ is recovered. Moreover, considering the low 
frequency approximation auj ^ 1 one finds that the reactive magnetic field is exactly the same as in the absence 
of field dissipation. Eliminating from the linearized equations ( ^ , ^ ) the gradient terms by operating twice with curl 
on the Navier-Stokes-equation and projecting on the z-axis we find 

{pdt - v^^)V^w = -aH,xtdtdzVH, + afioH^^^d^^w (80) 

^From the fact that iJj is harmonic. V^fe^ can be replaced by 

V^&,^MoM. , d.^'w. (81) 

1 — IfirCtOJ 

Using again the plane wave ansatz for the lateral dependence of w, the z-dependence is determined by 

{d'^ ~ k'){d! ~ f)w = , (82) 

with 

e = -— ^ (83) 

1 + K 

with the real part of q positive. The proper solution for w{z) remaining bounded for z ^ — cx) is hence 

W(z) = Ae*^^ + Be«^ (84) 

with the so far undetermined coefficients A and B. Using now aui <C 1 and expressing C, in terms of A and B by using 
the kinematic boundary condition the requirement of non-trivial solutions in A and B for the two stress boundary 
conditions results finally in (55). 
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FIG. 1. 

Sketch of the experimental setup. A ferrofluid of infinite lateral extension and infinite depth is exposed to a stationary 
magnetic field Hext directed normally to the surface. 
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FIG. 2. L 
(a) The decay rate to" {k,Hext) for free waves on the surface of a ferrofluid exposed to a normal magnetic Tiext- 
Solid lines (dots) denote the mesoscopic (hydrodynamic) theory. For the maximum investigated field intensity 
Hext = IhkA/m the relative deviation (b) does not exceed 4%. Parts (c,d) demonstrate how the wave excitation 
frequency uj' must be adapted to keep the the inviscid contribution ujQ{k, Hext) constant during the variation of the 
field intensity. For the evaluation of the curves we used Eqs.( p9| , |60[ ) and the material specifications for the ferrofluid 
APG J12 with an estimated magnetic relaxation time of r = 3.8 x lO^^s. 



